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ABSTRACT 


A  model  for  the  propeller-singing  phenomenon  considered  as  a  self- 
excited  oscillation  is  presented  to  Interpret  the  finding  of  a  recent 
experimental  work:  viz.,  that,  although  the  singing  frequency  roughly  obeys 
the  well-known  Strouhal  relation,  once  the  strong  singing  state  has  been 
established,  the  frequency  is  kept  constant  through  a  fairly  wide  range 
of  flow  velocity,  and  consequently  the  f requency-versus-veloci ty  diagram 
exhibits  step  and  Jump  characteristics.  The  model  presented  is  a  "closed 
loop"  composed  of  a  blade,  as  a  mechan ica l-v i brat  ion  system,  and  the  Karman 
vortex-shedding  mechanism;  the  blade  vibration  controls  the  shedding  mechan¬ 
ism,  and  the  hydrodynamic  reaction  of  shed  vortices  sustains  the  blade 
vibration.  The  control  imposed  by  the  blade  vibration  upon  the  vortex  shed¬ 
ding  actually  implies  the  synchronization  of  the  latter  with  the  former. 

The  model  which  simulates  the  vortex-shedding  mechanism  is  essentially  a 
simplified  mathematical  expression  for  the  disintegration  process  of  the 
vortex  sheets  shed  from  the  separation  points  into  the  rows  of  discrete 
vortices.  The  stability  criterion  derived  for  the  synchronized  run  of  the 
shedding  mechanism,  together  with  the  positive-work  criterion  imposed  upon 
the  phase  relation  between  the  blade  vibration  and  the  hydrodynamic  reaction 
of  the  shed  vortices,  gives  a  reasonable  interpretation  for  the  step  and 
jump  characteristics. 
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NOMENCLATURE 

shape  parameter  in  "influence"  functions  (Equation  [lO]) 
growth-rate  parameter  of  vortices 
form  constant  (Equation  [10]) 

shape  parameter  in  "influence"  functions  (Equation  [iO]) 
chord  iength 

integration  constant  (Equation  [20]) 
dimensionless  velocity  (Equation  [391) 
trai I ing-edge  thickness  of  body 

function  of  vortex- induced  velocity  (Equation  [391) 

function  defined  by  Equation  (31a) 

amplitude  of  hydrodynamic  reaction  force  on  blade 

function  (Equation  [26]) 

natural  frequency  of  blade 

shedding  frequency  of  Karman  vortices 

mode  shape  of  blade  vibration  (Equation  [46]) 


vii 


R-I059 


I 


H 

h 

k 


S,(x),  S2(x) 


S  -  (ud/U  -  2TTS, 


S'  -  (i)'d/U„  -  2nS' 

m  I 


'(.) 


%.) 


u. 


‘'o'  %(.) 


unit  step  function 

amplitude  in  blade-vibration  mode  shape 

characteristic  parameter  in  Hurwi t2-Routh  stability 
criterion  (Equation  [36]) 

order  of  iteration  process  for  criteria  of  stable 
synchronizat ion 

constant  >  I  or  2  (Equation  [39]) 

"influence"  function  (see  Equation  [27]) 

Strouhal  number  =  fj^ 
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6  phase  angle  between  flow  velocity  due  to  blade 

vibration  and  velocity  of  blade  vibration  at 
separation  point 


i  expression  given  by  Equation  (34) 


cp  angular  coordinate  along  blade  chord  (Equation  [40]) 

I 

X  expression  given  by  Equation  (26) 


¥  expression  given  by  Equation  (34) 


UD  angular  frequency  of  free  vortex  shedding 

angular  frequency  of  vortex  shedding  at  synchroni¬ 
zation  state  <■  natural  frequency  of  blade  vibration 


Subscripts 


I  refers  to  imaginary  part 

I  refers  to  lower  separation  point 


r 


refers  to  real  part 
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refers  to  upper  separation  part 
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refers  to  flow  velocity  induced  by  blade  vibration 
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INTRODUCTION 

The  propeller-singing  phenomenon  has  for  a  iong  time  been  understood 
as  a  forced  vibration  of  the  blade  due  to  the  hydrodynamic  reaction  force 
of  the  Karman  vortices  shed  from  the  traiiing  edge. 

Recentiy,  however,  an  experimental  study^  reveaied  a  particular 
feature  of  the  phenomenon:  though  the  relation  of  singing  frequency  versus 
flow  velocity  roughly  obeys  the  well-known  Strouhal  relation 

k  t  d 

the  frequency  in  the  strong  singing  state  Is  nearly  constant  through  a 
fairly  wide  range  of  flow  velocity  and  consequently  (as  is  shown  in  Figure  I) 
the  relation  of  singing  frequency  versus  flow  velocity  has  steps. 

Here 

f|^  shedding  frequency  of  the  Karman  vortices 
*  Strouhal  number 
■  flow  velocity 

d  ■  trai I ing-edge  thickness  of  body 
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Each  frequency  at  strong  singing  state  (characterized  by  a  step)  seems  to 

correspond  to  one  of  the  natural  frequencies  of  the  blade.  During  changes 

in  the  flow  velocity,  a  strong  singing  state  appears  on  each  step,  and 

successive  changes  cause  the  jump  phenomena  shown  by  the  dotted  lines  in 

the  figure.  This  kind  of  phenomenon  was  also  observed  in  the  experiment 

2 

on  suspension-bridge  oscillation  due  to  vortex  shedding. 

Arnold  et  al^  interpreted  this  particular  feature  of  the  strong¬ 
singing  phenomenon  as  a  special  kind  of  resonance  based  on  an  experimental 
hypothesis  that  the  shedding  frequency  of  the  Karman  vortices  becomes 
lower  with  increase  in  amplitude  of  blade  vibration.  However,  as  is 
suggested  by  Krivstov  and  Pernik,^  it  is  more  natural  to  regard  the 
singing  phenomenon  as  a  self-excited  oscillation  of  the  system  which  in¬ 
cludes  the  Karman  vort ices-shedding  mechanism.  In  this  paper,  a  model  for 
the  singing  propeller  is  presented  along  this  line. 
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SELF-EXCITED  OSCILLATION  LOOP  FOR  SINGING  PHENOMENON 

It  is  well  known  that  the  self-excited  oscillation  system  is  simu¬ 
lated  by  a  "closed  loop."  For  the  propeller-singing  phenomenon,  the  loop 
will  be  expressed  as  Is  shown  in  Figure  2.  The  "blade"  element  may  be 
regarded  essentially  as  a  mechanical-vibration  system.  The  other  element  — 
the  Karman  vortices-shedding  mechanism  —  Is,  on  the  other  hand,  a  self- 
excited  system  which  can  continue  to  shed  vortices  periodically  without  any 
periodic  stimulation  from  outside.  This  element  plays  the  more  important 
role  in  the  singing  phenomenon.  The  discrete  structure  of  the  Karman  vor¬ 
tices  indicates  that  this  self-excited  system  should  have  strong  non-linear 

characteristics.  The  important  aspect  of  such  a  non-linear  self-excited 

4 

system  Is  the  phenomenon  called  "synchronization"  or  "entrainment";  that 
is,  the  operation  of  the  non-linear  self-excited  system  is  often  synchron¬ 
ized  with  the  periodic  stimulation  from  outside  which  has  a  frequency  not 
so  different  from  the  natural  frequency  of  the  system.  In  the  present 
problem  the  above-mentioned  phenomenon  corresponds  to  the  synchronization 
of  shed  vortices  with  blade  vibration. 

The  blade  element  is,  as  stated  above,  essentially  a  mechanical- 
vibration  system  with  large  mass  and  stiffness,  and  definite  natural  fre¬ 
quencies  f.  .  The  natural  frequency  of  the  Karman  vortices-shedding 
®  U 

mechanism,  given  in  the  form  fj^  ■  depends  strongly  upon  the  flow 

velocity  .  in  the  range  of  In  which  fj^  differs  largely  from 

the  natural  frequency  of  the  blade  f|^  ,  the  vibration  amplitude  of  the 
blade  due  to  the  reaction  force  of  the  Karman  vortices  will  be  small,  and 
the  signal  from  the  blade  element  will  be  too  weak  to  synchronize  the 
Karman  vortices-shedding  mechanism.  In  this  condition,  the  loop  in  Figure  2 
is  open  between  these  two  elements,  and  the  system  is  in  the  state  of  a 
mere  forced  vibration  caused  by  the  reaction  of  the  shed  vortices.  This 
is  the  weak  singing  which  occurs  at  the  frequency  corresponding  to  the 
Strouhal  relation.  When  is  changed  and  f|^  approaches  f|^  ,  the 

amplitude  of  the  blade  vibration  becomes  large  enough  —  in  other  words,  the 
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signal  from  the  blade  element  becomes  strong  enough  —  to  synchronize  the 
shedding  of  the  vortices.  Thus,  the  loop  in  Figure  2  is  built  up  and  the 
system  enters  into  the  strong  active  state.  Of  course,  in  this  picture, 
it  is  necessary  that  the  loop  transfer  function  have  unstable  character. 

If  the  shedding  mechanism  of  the  vortices  can  be  synchronized  with  the 
vibration  of  the  blade  through  a  fairly  wide  range  of  velocity,  one  can 
see  that,  with  the  strong  singing  state  established,  frequency  remains 
constant  at  one  of  the  natural  frequencies  of  the  blade,  through  a  certain 
flow-velocity  range. 

The  strong  singing  state  of  the  present  model  will  occur  when  the 
Karman  vort Ices-shedding  mechanism  is  synchronized  with  the  blade  vibration 
and,  furthermore,  when  a  favorable  phase  relation  between  the  blade  vibra¬ 
tion  and  the  hydrodynamic  reaction  exists.  If  these  conditions  are  lost 
by  changing  the  flow  velocity,  the  system  will  Jump  from  the  strong  singing 
state  to  a  weak  singing  one  with  constant  Strouhal  number  and  will  continue 
in  the  latter  state  until  it  reaches  the  next  natural  frequency,  where 
another  cycle  of  strong  singing  will  appear. 
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THE  KARMAN  VORTICES  SHEDDING  MECHANISM  AND  ITS 
SYNCHRONIZATION  WITH  BLADE  VIBRATION 

To  provide  concrete  support  for  the  above  discussion,  a  suitable 
model  for  the  vort ices-shedding  mechanism  should  be  presented.  Von  Karman's 
Mork  notes  the  stability  of  the  vortex  rows,  but  it  does  not  explain  how 
these  rows  of  discrete  vortices  come  into  existence.  This,  together  with 
the  fact  that  we  lack  any  record  of  experimental  observations  of  the  flow 
near  the  trailing  edge  and  in  the  wake  of  an  oscillating  body,  proves  a 
stumbling  block  in  the  present  analysis.  A  mathematical  model  will  be 
presented  here  which  possesses  the  synchronization  mechanism  and  the  main 
features  of  the  above-mentioned  propeller-singing  phenomenon. 


PRESENTATION  OF  THE  MODEL 

The  vortex-shedding  process  may  be  dealt  with  by  solving  the  Navier- 
Stokes  equation  under  given  conditions.  The  analytical  approach,  however, 
seems  desperate,  and  the  numerical  attack  also  seems  hopeless  (except  for 
the  case  of  low  Reynolds  number)  even  if  a  computer  with  large  capacity 
Is  used.  Another  approach  may  be  to  treat  the  process  as  a  kind  of  Helmholtz 
instability  problem^  (using  the  method  adopted  by  Rosenhead^),  since  obser¬ 
vation  of  the  downstream  flow  of  a  body  with  a  blunt  trailing  shape  indicates 
that  the  vortex  sheets,  which  originate  from  the  vorticity  in  the  boundary 
layers  and  are  shed  from  the  separation  points,  roll  up  and  concentrate 
into  rows  of  discrete  vortices  a  rather  short  distance  away.  The  model 
presented  here  is  constructed  on  the  basis  of  this  description,  with  some 
speculative  mathematical  simplifications  which  are  adopted  in  order  to 
introduce  a  non-linear  oscillation  version  into  the  field  of  hydrodynamics 
and  to  obtain  a  wider  view  of  the  particular  feature  of  the  singing  phenome¬ 
non  described  in  the  previous  sections  of  this  report. 

As  noted  earlier,  the  vortex-shedding  mechanism  Is  a  self-excited 
system  and  itself  should  have  a  closed  loop,  which  plays  the  role  of  a 
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"minor  loop"  in  the  main  loop  shown  in  Figure  2.  On  the  basis  of  the  fore¬ 
going  description,  the  synthesis  of  the  loop  of  the  shedding  mechanism  may 
be  given  as  follows: 

(1)  The  existing  downstream  vortex  rows  exert  a  periodic  flow  distur¬ 
bance  upon  the  separation  points. 

(2)  This  flow  disturbance  at  the  separation  points  causes  a 
periodic  disturbance  in  the  strengths  of  the  vortices  which  are 
shed  from  the  separation  points  into  the  vortex  sheets. 

(3)  Thus  the  generated  non-uniformity  of  the  vorticity  distribu¬ 
tion  in  the  vortex  sheets  plays  the  role  of  the  embryo  of  the 
discrete  vortices;  or,  in  other  words,  this  non-uniformity  grows 
up  into  the  Karman  vortex  streets. 

The  block  diagram  of  the  thus  synthesized  loop  and  the  corresponding 
schematic  picture  of  the  flow  are  shown  in  Figure  3. 

If  the  blade  is  vibrating,  flow  disturbance  due  to  the  vibration 
may  be  superimposed  on  that  due  to  the  shed  vortices  at  the  separation 
points,  as  is  shown  in  the  block  diagram  of  Figure  3  by  a  dotted  line.  This 
flow  disturbance  due  to  blade  vibration  plays  the  role  of  the  synchroniza¬ 
tion  trigger  signal. 

MATHEMATICAL  EXPRESSION  FOR  THE  CONSTITUENT  ELEMENTS 

Assuming  two-dimensional  ideal-fluid  flow,  the  mathematical  expres¬ 
sion  of  the  constituent  elements  for  the  foregoing  model  may  be  given  as 
fol lows: 

(a)  The  Disturbance  Velocity  at  the  Separation  Points 
Induced  by  the  Downstream  Vortices 

The  velocities  induced  by  the  shed  vortices  at  the  separation  points 
are  affected  by  the  shape  of  the  body.  For  simplicity,  let  us  consider 
first  a  circular  cylinder  body  of  unit  radius.  The  coordinate  system  is 
shown  in  Figure  4a.  It  is  assumed  that  the  flow  separation  points  are 
located  at  (x  ■  0,  y  *  -  I)  and  that  the  vortices  shed  from  the  separation 


5 


R-1059 


points  are  flowed  away  along  th-  lines  y  ■  -  I  .  The  Induced  velocity 
U  (x)  of  an  Isolated  vortex  F  located  at  (x  -  x,  y  >  1)  and  of  Its 

corresponding  Image  (as  shown  In  Figure  4a)  upon  the  upper  separation  point 
(x  >  0,  y  *  I)  Is  given  as 

r 


(1) 


similarly,  the  effect  of  F^  at  (x  >  x,  y  >  -I)  and  of  Its  corresponding 
Image  upon  the  upper  separation  point  Is 


%<<■'>  ■  -  2S  ('  ■ 


(2) 


In  the  same  way,  the  effects  of  F^  and  F^  upon  the  lower  separation 
points  are  given  by 


“/„<»)  ■  2h('-7T?') 


(J) 


2tt 


(4) 


The  functions  .  Syjj(x)  ,  S^y(x)  ,  and  S^^(x)  are  Introduced, 


uu 


and  we  have 


S  (x) 
uu 


u  (x) 
uu 


S|jj(x) 


U|jj(x) 


(5) 


For  convenience,  these  are  called  "Influence  functions,"  since  they  represent 
the  effects  of  the  unit  strength  vortex  at  (x  -  x,  y  >  -  I)  on  the  Induced 
velocity  at  the  separation  points.  These  functions  are  graphically  ex¬ 
hibited  in  Figure  4b. 
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In  the  case  of  a  body  with  trailing  edge  of  paraboiic  shape  (Figure  Sa) 
given  by 

y*  “  I  -  (2/3)  X  (6) 


with  flow  separation  points  at  x=0,y--l,itls  assumed  that  the 
shed  vortices  are  flowed  away  along  the  lines 


which  are  tangential  to  the  body  surface  at  the  separation  points. 
The  conformal  transformation 


C 


(7) 


which  maps  the  flow  field  In  the  z-plane  Into  the  right  half  of  the 
C'plane  (Figure  5b),  provided  that  a  cut  along  the  x>axis  from  »“  to  4/3 
is  introduced,  determines  the  complex  velocity  potential  by  means  of  the 
"wall  effect."  The  influence  functions  S^^(x)  ,  S^jj(x)  ,  Sjjy(x)  , 
are  evaluated  and  exhibited  In  Figure  5c.  The  influence  functions  for 
both  circular  cylinder  body  and  parabolic  trai I ing-edge  body  exhibit  the 
following  general  features,  as  shown  by  Figures  4b  and  5c; 


<'»  V*  --htM 


m 


(2)  S  (x)  and  S..(x)  have  finite  value  at  x  ■  0,  while 

UU 

Syjj(x)  and  are  zero  at  x  -  0.  These  functions 

have  an  incubation  interval  before  a  large  rate  of 
increase  appears  with  increasing  x  ,  and  remain  almost 
constant  at  still  larger  x  . 
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(3)  For  large  values  of  x  , 


(9) 


By  utilization  of  these  general  characteristic  features,  the  influence 
function  can  be  approximately  expressed  by 

S  (x)  -  -  S,,(x)  «  -  Ae"® 

uu  XX'  ' 

Su^W  -  -  Ae""'^  +  Be‘^  H(b)  (iO) 

where  H(b)  is  the  unit  step  function  at  x  >  b  and  A  ,  B  ,  o'',  and 
B'  are  positive  constants. 

For  further  simplification,  the  motivation  of  which  will  become 
apparent  in  the  discussion  which  follows  in  a  later  section,  it  is  assumed 
that  B  ■  0  ,  ■  0,  and  a'  -  0  ;  hence 

"S/uW  ■  A  <"> 

This  simplification  affects  the  results  quantitatively,  but  presumably  not 
qual I  tat Ively . 


(b)  The  Strength  of  Shed  Vorticity  in  Terms  of 
Flow  Velocity  at  Separation  Points 

Let  the  velocity  distribution  in  the  boundary  layer,  the  velocity 
outside  the  boundary  layer,  and  the  thickness  of  the  boundary  layer  at  the 
upper  separation  point  be  designated  by  U  (y)  ,  ,  and  6  ,  respectively. 

The  vorticity  shed  per  unit  time  from  the  upper  separation  point  Into  the 
upper  vortex  sheet  is  given  by 


/**  d  U  (y)  1 

'su  "  X  ^  dy  -  ■?  Uf. 


dy 


(12) 
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Similarly,  for  the  lower  separation  point, 

'■si  ■  1“!  <'2*) 

where  is  the  velocity  outside  the  boundary  layer  at  the  lower  separa¬ 

tion  point. 

If  it  is  assumed  that  the  velocity  at  the  upper  separation  point, 

U  ,  Is  composed  of  the  mean  velocity  U  and  a  small  perturbation  velocity 

^  i  U)  t  ^  ^ 

of  periodic  nature,  »  then  the  vorticity  shedding  rate  from  the 

upper  separation  point  is  given  by 

r  fu®  +  2U  u  (13) 

su  2  \  ou  ou  ou  I 

and  similarly  from  the  lower  separation  point,  by 


(c)  Concentrating  Process,  Vortex  Sheets  into 

Discrete  Vortex  Rows 

Rosenhead's  treatment^  for  the  rolling-up  phenomenon  of  the  vortex 
sheet  In  the  Helmholtz  instabi 1 1 ty  problem  indicates  that  the  local  distur¬ 
bance  in  the  vorticity  density  in  the  vortex  sheet  causes  a  geometrical 
deformation  of  the  sheet.  This  deformation  gives  rise  to  the  migration 
of  the  vortices  in  a  manner  which  promotes  the  growth  of  vorticity  density; 
thus  these  two  processes  cooperate  in  causing  an  accelerating  concentration 
of  vorticity.  In  the  present  treatment,  a  simple  mathematical  expression, 
based  on  the  above-mentioned  feature  of  the  Helmholtz  instability  problem, 
is  given  for  the  disintegration  process  involving  the  shedding  of  the  vortex 
sheets  from  the  separation  points  into  discrete  vortex  rows. 

It  is  assumed  that  the  process  follows  the  pattern  outlined  below 
and  shown  in  Figure  6. 

(I)  The  vortex  sheets  shed  from  the  separation  points  are  flowed 
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away  along  the  designated  path. 

(t 

(2)  The  vorticity  densities  in  the  sheets  have  sinusoidal  distur¬ 
bances  due  to  the  periodic  fluctuations  in  the  vortlcity- 
shedding  rate  given  by  Equations  (13)  and  (|4), 

(3)  The  amplitudes  of  the  disturbances  Increase  with  distance 
X  from  the  separation  points. 

According  to  this  picture  the  vorticity  distributions  in  the  upper  and  the 
lower  sheets  may  be  expressed  as 

lu)tt  -  (x/UjJ 

r  (x.t)  -  r  +  Y  (x,t)  -  r  +  y  (x)e 
u'  ’  ou  u'  *  '  ou  ou 

and  (15) 

Itult  -  (x/U  )1 

respectively,  where  is  the  constant  velocity  in  the  x-direction  by 

which  the  vorticity  Is  flowed  away,  T  and  f  ,  are  the  mean  vorticity 
densities  in  the  sheets,  and  YQy(x)  and  V^^(x)  are  the  amplitudes  of 
the  disturbance  terms  in  the  vorticity  densities.  On  the  basis  of  the 
continuity  equation 

r  At  -  r  u  At 

su  u  w 


where  At  is  a  specified  Interval  of  time,  together  with  Equations  (13) 
and  (14),  It  Is  apparent  that 


"ou 

1 

*  2 

o-^ 

c 

o 

U  u 
ou  ou 

“  U 

w 

(16) 

and 

"oi 

-1 

“  2 

‘^oi* 

•'w 

-U  ,u  , 
ol  ol 

“  U 

w 

(17) 

It  will  be  assumed  that  the  process  which  embraces  the  passing  of 
the  vortices  from  the  embryonic  stage  to  complete  disintegration  into 
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discrete  vortices  Is  completed  at  the  stage  where  the  amplitude  of  the 
perturbation  vorticity  density  becomes  equal  to  the  mean  vorticity  strength 
of  the  sheet. 

The  simplest  mathematical  expression  which  satisfies  the  above 
requirement  and  fixes  the  growing  rate  with  x  of  the  disturbance  ampli¬ 
tudes  I  Y^yCx)!  and  |  Yq£(x)|  is  assumed  to  be 


l^o' 


Yo(x)| 

-  (1 


'"o' 


I  "o' 


(18) 


where  the  constant  a  is  the  growth-rate  parameter.  The  absolute  value 
has  been  introduced  so  that  the  above  relation  will  hold  true  for  either 
upper  or  lower  vortex  sheet.  It  Is  obvious  that  this  relation  ensures  the 
final  value  of  |  Yq(x) |  as  equal  to  |r^|  . 

The  solution  is  given  by 


I  V^(«)l  ,axtc 
|rj  ■ 


where  the  constant  c  is 


c 


IvqM^oI 
I  -  I  v„(o)/rJ 


(19) 


(20) 


The  value  of  YQ(x)/|r^|  Is  graphically  exhibited  (see  dotted  lines)  In 
Figure  7,  in  terms  of  ax  +  c  .  If  for  simplicity  the  value  of 
|Yo(x)|/|r^|  is  replaced  by  a  unit  step  function  at  ax  +  c  ■  0  ,  then 
It  can  be  seen  (from  Figure  7)  that  a  suitable  approximation  of  the 
function  lYQ(x)|/|r^|  is  obtained.  Hence 


Iy^(x)| 


H 


(21) 
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where  H(x)  denotes  a  unit  step  function  at  x  .  Applying  Equation  (21) 
to  Equations  (15),  (16),  and  (17),  the  vorticity  distribution  on  the  upper 
and  lower  vortex  sheets  will  be  given  by 


and 


U  .  -f  u  .e 
OJt  ol 


iu)t 


are  the  velocities  at  upper  and  lower  separation  points  respectively. 


FREE  SHEDDING  OF  THE  KARMAN  VORTICES 

Now  it  is  possible  to  cot.  ^ose  a  mathematical  equation  for  the  closed 
loop  of  the  vort iceS'Shedding  mechanism.  The  closed-loop  equations  with 
respect  to  the  upper  and  the  lower  separation  points  are 

-  Uco  ■  jT  r^(>‘.t)  +  s^^(x)  r^(x,t)j  dx  (23) 

where  denotes  the  uniform  flow  velocity  at  infinity  and  the  other 
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notations  are  as  given  In  Equations  (9)  and  (lO). 

Here,  the  following  assumptions  are  Introduced: 

(a)  Y^^(x)  -  Y^/x)  (24) 

Therefore,  from  Equations  (16)  and  (I7)t 


u  -  -  u  , 
ou 

As  can  be  seen  In  Figure  6,  this  assumption  Is  imposed  by  the  well-known 
arrangement  of  vortices  typical  to  the  Karman  vortex  street. 

(b)  U  -  U  ,  -  -  2U  (25) 

'  ou  oX  “>  W 


If  the  body  Is  flat  in  the  direction  of  the  main  flow,  and  if  the  wake 
region  near  the  trailing  edge  is  regarded  as  dead  water,  the  relations  of 
Equation  (25)  should  give  a  good  approximation. 

Then  Equation  (23),  together  with  Equation  (22),  yields  the  time- 
dependent  part  of  the  shedding  mechanism: 


u  e 
ou 


Iwt 


u 

ol 


lu)t 


where 


■f:\ 


s  (x)  — sa-  H 
““  (u  I 

ou 


luu(t-2x/U  )  . 
e  ®  dx 


ol 


'  I VI 

i«>(t-2x/U„) 


f(x)  »  jtn  and  x 


2u 


ou.Jt 


(26) 
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If  the  symmetrica)  relations,  Equations  (6)  and  (24),  are  taken  Into 
account,  the  two  equations  In  Equation  (26)  become  Identical  and  can  be 


expressed  as 

u 

o 

" 

rh- 

j].-"-  dx 

(27) 

where 

u  - 

o 

U  ■  -  U  J 

ou  o« 

S,(x) 

- 

-  -  S^^(x)  -  llm  j 

Ss(x) 

■ 

..s,„(x)  .  I!»  ( 

Qf-O  ‘ 

-Ae"“  *H(b)j 

S  - 

2(ii/U 

00 

If  the  representative 

thickness  of 

the  body  Is  taken  as 

2  (the  value 

taken 

as  the  separatlon«polnt  thickness) 

In  the  evaluation  of 

(x)  and 

S,(x)  . 

then 

S 

-  2ttS 

(28) 

where  Is  the  Strouhal  number. 

It  should  be  noted  that,  due  to  the  symmetry  relations  Introduced 
by  Equations  (8),  (24),  and  (25),  the  system  of  simultaneous  equations 
given  by  Equation  (23)  or  (26)  for  the  upper  and  lower  separation  points 
Is  reduced  to  a  single  equation  (Equation  [27]). 

Upon  integration.  Equation  (27)  yields 
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or 


J. 

a 


f 


The  real  and  Imaginary  parts  of  the  above  expressions  are  sufficient  to 
determine  the  unknown  S  and  -rp  . 


SYNCHRONIZATION  OF  SHEDDING  MECHANISM  WITH 
BLADE  VIBRATION  -  STABILITY  CRITERIA 


If  the  blade  vibration  Induces  periodic  flow  disturbances  u  e^**’  ^ 
luj't 

and  u^^e  at  the  upper  and  the  lower  separation  points  respectively, 

these  may  act  as  a  synchronization  signal,  as  Is  shown  In  the  block  diagram 
of  Figure  3.  Assuming  the  symmetrical  relation 


vu 


u  ) 

V 


(30) 


the  loop  equation  for  the  synchronized  condition  is  obtained  as 


OD 


(31) 


where 


s'  -  u)' -2- 


and  a  is  the  phase  relation  between  u^  and  u^  . 

It  should  be  noted  that  in  the  above  equation  the  periodic  disturbance 
induced  by  the  blade  vibration  is  assumed  to  be  known,  whereas  the  total 
perturbation  velocity  at  the  separation  point,  u^,  and  the  phase  angle  or 
are  both  unknown  and  will  be  determined  by  the  solution  of  Equation  (31) • 

It  should  be  pointed  out,  also,  that  Equation  (31)  describes  the  neutral 
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condition  for  the  synchronization,  and  that  there  is  thus  a  possibility 
that  Its  solution  will  represent  the  unstable  synchronization  state,  where 
small  deviations  from  the  neutral  condition  may  grow  unboundedly  with  In¬ 
creasing  time.  To  exclude  this  possibility,  a  stability  analysis  for  the 
synchronization  must  be  performed.  The  powerful  method  introduced  by 
Van  del  Pol  for  stability  analysis^  cannot  be  utilized  in  the  present  case, 
since  Equation  (31)  Is  an  algebraic  equation  stating  simply  the  neutral 
condition;  hence  the  timewise  growth  and  decay  of  any  deviation  from  this 
condition  Is  not  taken  into  account.  In  this  paper,  the  following  approxi¬ 
mate  method  will  be  used. 

An  Iteration  scheme  is  developed  based  on  the  fact  that  the  n^^ 
Iterative  values  of  u^  and  a  ,  written  In  the  form 


u  >  u  -f  (Au  ) 
o  oo  on 


O'  -  Of  +  (AOf) 
o  n 


can  be  obtained  from  their  (n-l)th  values. 


and 


u  -  u  +  (Au  )  , 

o  oo  on- 1 


»  ■  “o  * 


(where  u  and  a  are  the  values  at  the  neutral  condition  and  Au  and 
'  oo  o  o 

Aof  are  small  deviations  from  these  values)  by  means  of  Equation  (31).  Let 


Then  Equation  (31)  can  be  written  as 


u 


-2 


lo 

e 


(31a) 
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Substitution  of  the  Iterative  values  on  the  left-hand  side  and  the 
(n-))th  values  on  the  right-hand  side  of  this  equation  yields 


u  + 
_20 _ 


(Au  ) 

o  n 


If  only  the  zero-  and  first-order  terms  are  kept  in  this  iteration  scheme, 
the  following  approximate  equation  results: 


From  Equation  (31a),  at  the  neutral  condition. 


u  lor 
oo  o 


Substitution  of  Equation  (33)  into  Equation  (32)  yields 


(33) 


—  *  I 


o  n* 


oo 


(Aor) 


n-l 


Subtracting 


(Au  )  , 

o  n- 1 
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from  both  sides  of  the  foregoing  equation  results  in 


where 


(34) 


i 


u 

oo  °° 
“  oo 


Since  (Au  )  >  (Au  )  ,  and  (Ao')  -  (Acv)  ,  are  the  changes  in  Au  and 

'  o  n  '  o  n-1  '  'n  '  'n-1  ^  o 

Aof  per  unit  operation  from  n-1  to  n  ,  or  the  rate  of  change  of  Au^ 
and  tot  ,  Equation  (34)  can  be  written,  for  the  case  of  very  slow  rate  of 
change,  as 


J_ 


d(Au^) 

dn 


+  I 


*^oo  d(Acr) 
dn 

OD 


(35) 


which  represents  a  system  of  first-order  homogeneous  differential  equations 
with  unknown  Au^  and  Aar  .  On  assuming  solutions  of  the  form 


Au 


u  e 
o 


kn 


and 


oo 


(4a) 


X  e 
o 


kn 


the  following  characteristic  equation  is  obtained 


(36) 
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where  r  and  I  denote  real  and  imaginary  parts,  respectively.  The 
growth  or  decay  of  Au^  and  AQ^  is  discriminated  by  applying  the  Hurwitz- 
Routh  criteria.  For  a  stable  solution,  i.e.  decay  of  the  deviations  from 
the  neutral  state,  the  conditions 


(a) 

-(f^  +  «,)  >  0 

(b) 

»r<p,  -  >  0 

(37) 


must  be  met . 

The  previously  raised  question  of  the  stability  at  the  synchronization 
stage  can  now  be  tackled  on  the  assumption  that  there  is  a  very  slow  time- 
wise  change  of  the  deviations  Au  and  Acr  .  In  such  case  the  problem  of 
stability  at  synchronization  is  equivalent  to  the  problem  of  the  convergence 
of  the  foregoing  Iterative  process.  Then  Equation  (37)  nay  be  regarded  as 
the  stability  criterion  for  synchronization. 

Upon  integration,  as  In  Equation  (29),  Equation  (31)  yields 


or 


“o  lor 

u:*  * 


if  -'S''’]."' 

»>1,(A) 


for 


(38) 


and  the  corresponding  stability  criterion.  Equation  (37),  with  the  use  of 
(38),  is 

(I)  I  -  Eq  cos  [s'  ^  f  (^)]+  D  cos  >  0 


19 


R-1059 


(n) 


where 


(39) 
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SYNCHRONIZATION  SIGNAL  FROM  THE  BLADE  VIBRATION  AND 
HYDRODYNAMIC  REACTION  OF  THE  SHED  VORTICES 

In  order  to  complete  the  closed  loop  shown  in  Figure  2,  the  synchron¬ 
ization  signal  due  to  the  blade  vibration  and  the  hydrodynamic  reaction  of 
the  shed  vortices  upon  the  blade  should  be  evaluated.  In  the  mathematical 
treatments  of  those  two  quantities,  the  blade  is  assumed  to  be  a  flat  plate, 
and  the  thin-airfoil  method  is  used.  The  corresponding  coordinate  system 
Is  shown  in  Figure  8. 


SYNCHRONIZATION  SIGNAL  DUE  TO  BLADE  VIBRATION 


The  flow  in  the  system  under  consideration  is  composed  of  two  parts. 
One  is  "no  circulation"  flow  due  to  the  blade  vibration  and  the  other  is 
the  flow  with  circulation  due  to  the  shed  vortices.  "No  circulation"  means 
that  the  total  circulation  around  the  body  is  zero.  In  Equation  (31)  or 
(38),  that  flow  velocity  at  the  separation  points  which  is  due  to  the  "no 
circulation"  part  is  u^  ,  and  that  due  to  the  total  flow  Is  u^  .  Deriva¬ 
tion  of  u^  for  the  assumed  mode  of  blade  vibration  follows. 

Consider  the  bound-vortex  distribution 


Y(,(x,t) 


luj't 

e 


2U_  A 


cos  <P  _  1  . 
s  In  <p  2  1 


cos  2cp 
s In  9 

n-2 


A 


n 


sin  ntp 


] 


itu't 

e 


r 

where  x  -  -  — (I  +  cos  tp)  and  C  is  the  chord  length.  Clearly,  this 
bound-vortex  distribution  satisfies  the  "no  circulation"  condition 


(40) 


(41) 


The  strength  of  the  bound-vortex  distribution  near  the  trailing  edge 
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9  ■  tt  -  e  ,  where  n  >  >  e  >  0  ,  is  given  by 

.  n  .  e)  ~  -  2UJA^*iA,)  i  (M) 

Therefore,  if  <p  ■  n  -  e  denotes  the  coordinate  of  the  flow  separation 
points, 

u  -  -  u  .  -  -  U„(A  +  4  A,  )  7  (^3) 

VU  v£  O  2  1  € 

This  satisfies  the  symmetrical  relation  Equation  (30),  and  u^  can  be 
written  as 

“v  ■  -  ".<*0  W  *.  >  e 

For  a  particular  vibration  mode,  A  and  A^  are  determined  as  follows: 

The  Induced  velocity  of  the  bound  vortex  in  the  y-direction  is 

v,(xt) 


If  the  vibrational  velocity  of  the  blade  in  the  y-dIrection  is  given 

Vjj(x,t)  =  U„g(x)  e""  (46) 

then  the  boundary  condition  on  the  blade  requires  that 

-  (^*"1;)  ''b  ■  “ 

where  y^^  is  the  vibrational  displacement  of  the  blade.  The  above 
boundary  condition  can  be  written  as 

V, (X, t)  -L(x,t)  +  Vj,(x,t)  dtl  -  0 

^  n  ^ 


\o^^^  lu>'t 


A  +  ZA  cos  ncp 
on 


,1  '‘"'t 

I 


(45) 
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and  therefore  and  Aj^  are  determined  by 

-  A^  +  2_^  A^  cos  ncp  =  g(x)  -  i  •—  g'(x)  (47) 

n=1 

HYDRODYNAMIC  REACTION  FORCE  ON  THE  BLADE 
DUE  TO  THE  SHED  VORTICES 

To  ascertain  the  existence  of  self-exciting  singing,  the  vibratory 
force  exerted  on  the  blade,  due  to  the  bound  vortex  and  to  shed  vortices, 
has  to  be  evaluated  and  its  phase  compared  to  the  vibration  velocity  of 
the  blade. 

The  hydrodynamic  force  exerted  by  the  "non-c i rculatory"  component 
of  the  bound-vortex  distribution  corresponding  to  any  blade  vibrational 
mode  has  a  virtual-mass  term  only  and  no  damping  component.  On  the  other 
hand,  the  reaction  force  of  the  shed  vortices  has  a  damping,  as  well  as  a 
virtual -mass,  component.  Therefore,  for  the  discussion  on  the  existence 
of  self-exciting  vibration,  only  the  latter  hydrodynamic  force  should  be 
taken  into  account. 

The  reaction  force  due  to  the  shed  vortices  is  evaluated  by  the 

/  /  8 

method  used  by  Karman  and  Sears  in  non-stat ionary  wing  theory.  As  is 

shown  in  Figure  8,  the  shed  vortices  are  assumed  to  be  flowed  away  along 
the  x-axis.  The  shed  vortices  have  their  counter  vortices  on  the  blades; 
these  are  known  as  bound  vortices.  Thus  the  flow  system  due  to  the  shed 
vortices  is  composed  of  the  vortex  pairs.  Then  the  momentum  In  the  y-dlrec- 
tion  due  to  the  i^^  vortex  pair  is  given  by 

I  *  pr.(x^j-x^,)  (48) 


where  p  ,  f.  ,  and  x,  .  denote  the  density  of  fluid,  the  circulation  of 
’  I  bi 

the  shed  vortex  located  at  location  of  the  counter  vortex, 

respectively.  The  force  upon  the  blade  due  to  the  particular  pair  is 


F 


dt 


’'bi) 


(49) 
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Assuming  for  simplicity  that  the  counter  vortex  Is  concentrated  at  the 
trailing  edge  (x  -  0) ,  Equation  (49)  becomes 


^ 


(50) 


In  the  present  problem,  the  sum  of  the  vorticity  distributions  from  the 
upper  and  the  lower  separation  points  yields 


r,  -  2e‘XM[i  ‘  (51) 


where  the  time  origin  is  the  same  as  in  the  case  of  Equation  (31)  •  The 
total  force  acting  upon  the  blade  is  then  given  by 


F,  --P 


a> 


2pU„u)' 


where 


r,xdx  - 


"  ■  i  'ft) 


n. 


^  1  (cr  +  p  -  tan  3  +  -j) 


(52) 


P  ■  tan  ^wS ' 


If  the  effective  thickness  of  the  blade  at  the  separation  point  is  taken 
as  d  , 

u)  d 


S  ■ 


U_ 


and  hence 


^ 


m 


+  vr  e 


1  (or  +  P  -  tan  P  +  -j) 


(53) 
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EXISTENCE  OF  THE  SELF-EXCITED  SINGING  STATE 


The  various  constituent  elements  of  the  closed  loop  comprising  the 
propeller-blade  vibration  and  the  Karman  vort Ices-shedd Ing  mechanisms  having 
been  determined,  it  Is  now  possible  to  present  a  mathematical  expression  for 
the  existence  of  the  self-excited  singing  state.  In  doing  so,  however,  It 
seems  useful  to  summarize  the  previously  obtained  results. 

(1)  Blade  Vibration  to  Synchronizat ional  Signal 

When  the  vibrational  velocity  of  a  point  on  the  blade  Is  given  by 

Vj^(x.t)  ■  U^g(x)  e'^  ‘  (46) 


the  synchronization  signal  Is  expressed  In  the  form  of 


•“>1  „  ^  1  A  \  •  .•‘U  t 

Uy  « 


(44) 


where 


-  A. 


COS  nq)  -  g(x)  -  I  g'(x)  ,  x  -  -  |  (1  +  cos  cp)  (4?) 


and  (p  ■>  TT-e(TT>>e>0)  is  the  location  of  the  separation  points. 


(2)  Synchronized  Operation  of  the  Shedding  Mechanism  with  u^e**^  * 

The  governing  equation  for  the  synchronization,  after  use  is  made 
of  the  simplified  expression  Equation  (11),  is  given  by 


%  la 


for 
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or 


Furthermore,  the  stability  of  the  synchronization  state  requires  that 
(I)  1  -  Eq  cos  fs ^  “  f 


(^)] 


+  0  cos  a  >  0 
o 


and 


(II) 


(39) 


Equations  (38)  and  (39)  give  u^  and  o'  in  the  stable  synchronized 
condition  for  given  synchronization  signals  u^  and  S'. 

(3)  Hydrodynamic  Force  Reaction  Due  to  the  Shed  Vortices 

When  u  and  a  are  given,  the  reaction  force  on  the  blade  Is 


-#p 


i  (or  +  P  -  tan  P 


+  ,) 


(53) 


where  the  origin  of  time  i.'  taken  so  that  u^  has  a  real  positive  value. 

If  the  above  results  are  used,  the  criterion  for  the  positive  work 
of  the  reaction  force  upon  the  blade  vibration  can  be  obtained;  this  is 
composed  of  the  criteria  for  the  existence  of  the  self-excited  singing 
state  together  with  that  of  stable  synchronization,  Equation  (39),  I  or  ii. 

In  the  first  place.  It  may  be  assumed  that  the  reaction  force  is  concentrated 
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near  the  trailing  edge  of  the  blade.  As  was  stated  earlier  (p.  23)  the 
"no  circulation"  component  of  the  bound  vortices  which  Is  caused  by  the 
blade  vibration  exerts  no  damping  effect.  Therefore,  If  the  absolute  value 
of  the  phase  difference  between  the  reaction  force  ^  and  the  vi¬ 

bration  velocity  of  the  trailing  edge  V|^(o,t)  -  U^g(o)e*'‘’  ^  Is  less  than 
,  the  reaction  force  should  exert  positive  work  upon  the  blade  vibration. 
As  Equation  (53)  has  Its  phase  origin  at  u^  ,  the  cr I  ter  ion  for  positive 
work  Is  given  by 


-TT<o+P-tanP+e<0 


(54) 


where,  from  Equations  (46)  and  (44), 


e  - 


tU)'t 

u  e 

■'"9  Vjj(9-"-e,t) 


-  arg 


-Ky  j  *.)j 

g(9-n-e) 
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NUMERICAL  CALCULATIONS  AND  DISCUSSION 

/  / 

FREE  SHEDDING  OF  KARMAN  VORTICES 

For  several  values  of  A  ,  b  ,  and  a  ,  the  unknown  Strouhal  number 
and  dimensionless  velocity  ratio  are  determined  by  Equation  (29) 

for  the  case  of  free  shedding  of  Karman  vortices.  Both  these  unknowns  are 
determined  by  solving  the  real  and  imaginary  parts  of  the  above  equation 
simultaneously.  Due  to  the  complicated  form  of  the  function  f(2u^/U„) 
the  equations  are  solved  graphically  and  the  corresponding  results  are 
sumtnari::ed  in  the  following  table. 


'•  a 

A 

u  /U 

O  0® 

2 

l/2n 

0.202 

0.141 

2 

0.8/2TI 

0.197 

0.127 

1 

0.6/2n 

0.175 

0.145 

2 

0.6/2TT 

0.191 

0.110 

3 

0.6/2TT 

0.200 

0.080 

2 

0.5/2TT 

0.189 

0.099 

The  parameters  A  and  b  ,  which  are  shape  parameters,  can  be  sel¬ 
ected  according  to  the  particular  body  shape.  The  parameter  a  ,  however, 
which  is  related  to  the  rate  of  growth  of  the  disturbance  in  the  vortex 
sheet,  can  only  be  selected  by  trial.  The  selected  values  of  A  and  b 
seem  to  be  reasonable  ones,  since  (A  ■  •  b  >  2)  and  (A  ^  •  b  ■  2) 

describe  approximately  the  circular  cylinder  body  and  parabolic  trailing- 
edge  body.  The  corresponding  value  of  Strouhal  number  (around  0.2)  seems 
very  close  to  the  usual  experimental  results.  The  values  of  seem 

reasonable  also.  This  is  an  indication  that,  in  spite  of  all  the  assumptions 
made,  the  present  mathematical  model  for  the  free-shedding  mechanism  is 
promising. 
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SELF-EXCITEO  SINGING 

A  series  of  numerical  calculations  was  made  to  Illustrate  the 
existence  of  the  self-excited  singing  and  to  clarify  Its  mechanism.  As 
Is  seen  In  Figure  9a,  the  following  vibrational  nndes  of  the  blade  have 
been  considered: 

(a)  1st  mode:  spanwise  bending  mode,  with  g(x)  >  h  ■  constant 

(b)  2nd  mode:  torsional  mode,  with  g(x)  -  -  h  cos  <P 

(c)  3rd  mode:  chordwise  bending  mode,  wlth'g(x)  ■  h  cos  2<p 

where  the  transformation  from  Cartesian  to  angular  coordinates  Is  given 
by  X  ■  -  ‘^(l+cos  <p)  ,  C  -  chord  length  of  blade.  The  ratio  of  the  separa¬ 
tion-point  thickness  d  to  the  chord  length  C  Is  0.1. 

The  numerical  work  has  been  performed  along  the  line  Indicated  In  the 
section  on  existence  of  the  self-excited  singing  state  (p.  29).  In  the 
first  place,  the  real  and  Imaginary  parts  of  Equation  (38)  are  solved 
simultaneously  with  respect  to  and  or  for  given  u^/Ugo  and 

It  may  be  useful,  here,  to  note  the  following  points: 

(1)  ,  the  trigger  signal  for  the  synchronization  of  the 
vortex  shedding  with  the  blade  vibration,  is  proportional 
to  the  amplitude  of  the  blade  vibration  and  Is  nothing  but 
the  amplitude  of  the  blade  vibration  expressed  in  hydro- 
dynamical  terms. 

(2)  -  tt>*d/U„  Is  the  non-dimenslonal  expression  of  the 
natural  frequency  of  the  blade  vibration  ou'  (■  the 
strong  singing  frequency)  In  terms  of  the  blade  thick¬ 
ness  at  the  separation  points,  d  ,  and  the  flow  velocity 
at  Infinity,  .  The  expression  S Vs  will  be  used  In 
place  of  s',  where  S  -  tt>d/U„  ,  and  ou  Is  the  free- 
shedding  angular  frequency  of  the  vortices.  Since 

s'/S  -  w*/w  ,  sVs  should  be  named  "tuning  factor." 

Meanwhile,  since  S  -  2ttS^  (see  Equation  f28] ;  S^  is 
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the  Strouha)  number)  is  constant  for  a  given  blade  shape 
independent  of  the  flow  veiocity,  and  since,  furthermore, 
u)'  and  d  are  also  constant  for  the  given  blade,  S  Vs  ■ 

((u^d/2nS^)  (1/U^)  has  the  sense  of  a  non-dimensional  ex¬ 
pression  of  the  reciprocal  of  the  flow  velocity  . 

If  the  obtained  values  of  and  a  are  applied  to  Equations  (39)'I 

and  (39)-IIj  the  limits  are  determined  for  the  stability  of  the  synchroniza¬ 
tion  of  the  vortex  shedding  with  the  blade  vibrat ion,  under  the  given  blade- 
vibration  amplitude  and  the  given  tuning  factor  sVs  .  In  a 

similar  way.  Equation  (5^)  can  be  used  to  check  and  determine  whether  or 
not  the  hydrodynamic  reaction  of  the  shed  vortices  exerts  positive  work 
upon  the  blade  vibration. 

In  the  actual  calculation,  owing  to  the  complicated  functional  form 

of  ^  difficult  to  solve  Equation  (38)  with  respect  to 

Uq/U„  and  a  for  given  and  S'/S  ,  and,  therefore,  an  auxiliary 

procedure  Is  introduced:  Equation  (38)  is  solved  with  respect  to  u^/U^ 
and  ot  for  assumed  u  /U^  and  sVs  .  This  procedure  gives,  for  each 
value  of  S Vs  ,  the  auxiliary  diagrams  illustrated  in  Figures  9b  and  9c, 
on  which  are  also  indicated  the  ranges  where  the  criteria  of  Equations  (39) 
and  (34)  are  satisfied.  From  these  auxiliary  diagrams,  the  final  charts 
shown  in  Figures  10a,  b,  and  c  are  derived.  Figures  10a,  b,  and  c  exhibit 
the  two  kinds  of  limit  boundaries  corresponding  to  the  stable-synchroniza¬ 
tion  criterion.  Equation  (39)-II,  and  the  positive-work  criterion.  Equa¬ 
tion  (54),  on  a  versus  S'/S  plane.  Since  the  condition  imposed 

by  Equation  (39)"II  Js  always  more  severe  than  that  imposed  by  Equation  (39) “I » 
the  limit  boundary  due  to  the  latter  is  omitted. 

The  values  used  in  this  calculation  for  the  parameters  are  A  ■  , 

b  ■  2,  and  a  »  I,  2,  3-  The  first  two  are  approximating  values  for  a 
parabolic  shape  of  the  trailing  edge  of  the  body.  On  the  other  hand,  the 
given  values  of  a  are  mere  trial.  However,  although  the  choice  of  a 
causes  some  quantitative  effect,  as  seen  in  Figures  lOa,  b,  and  c,  the 
qualitative  nature  does  not  seem  to  be  altered. 

It  should  be  noted  that  Figures  lOa,  b,  and  c  indicate  whether  or  not 
the  assumed  value  of  satisfies  the  two  criteria,  but  give  no 
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Inforinat  Ion  with  regard  to  the  determination  of  .  This  is  the 

problem  of  the  limit-cycle  amplitude  of  the  blade  vibration,  which  should 
be  resolved  by  considering  the  balance  between  the  energy  input  due  to 
the  hydrodynamic  reaction  force  from  the  shed  vortices,  F|^e*^  and  the 
energy  dissipation  due  to  the  damping  capacity  of  the  blade.  In  regard 
to  this  point,  it  should  be  repeated  that  the  non-circulation  flow  compo¬ 
nent  due  to  the  blade  vibration  exerts  no  damping  effect,  and,  therefore, 
only  the  mechanical  damping  should  be  considered  as  the  damping  capacity 
of  the  blade.  The  existence  of  the  limi'^  cycle  should  be  proved  through 
a  strict  mathematical  consideration,  but  it  may  be  anticipated  by  the 
fact  that  the  magnitude  of  the  reaction  force  is  not  affected  by  the 

amplitude  of  the  blade  vibration  as  much  as  the  mechanical  damping  force 
Is  affected.  With  increasing  amplitude  of  blade  vibration,  the  reaction 
force  remains  almost  unchanged,  whereas  the  damping  force  increases  so 
that  an  energy  balance  due  to  these  forces  will  eventually  appear.  If  the 
mechanical  damping  of  the  blade  is  known,  one  may  obtain  the  locus  of  the 
limit-cycle  amplitude  on  a  u^/U^  versus  S Vs  plane  (that  is,  an  ampli¬ 
tude  versus  reciprocal-of-the-f low-velocity  plane).  If  the  locus  exists 
in  the  region  bounded  by  the  two  kinds  of  criterion  limits  (the  region 
indicated  by  mark -•■■■  in  Figures  lOa,  b,  and  c)  through  a  certain 

range  of  S ^/S  ,  then  the  self-excited  strong  singing  which  manifests  It¬ 
self  as  a  step  on  the  frequency-versus-veloci ty  diagram  may  appear  through 
the  corresponding  range  of  the  flow  velocity.  When  the  limit-cycle  locus 
cuts  the  stable-synchronization  limit,  the  synchronization  of  shed  vortices 
with  blade  vibration  may  cease.  In  other  words,  the  closed  loop  shown  In 
Figure  2  is  opened,  and  consequently  there  should  be  a  jump  from  the  self- 
excited  strong  singing  state  to  the  weak  one  with  the  Strouhal  frequency. 
While,  as  seen  in  Figures  lOa,  b,  and  c,  the  stable  synchronization  limit 
seems  to  occur  always  as  a  lower  boundary  below  which  the  stability  is 
lost,  the  positive-work  limit  manifests  itself  either  as  an  upper  or  a 

lower  boundary,  above  or  below  which,  respectively,  the  favorable  phase 

iou^t 

relation  between  the  blade  vibration  and  the  hydrodynamic  reaction  F|^e 
vanishes.  When  this  limit  makes  an  upper  boundary,  the  limit-cycle  locus 
does  not  go  up  across  the  limit,  because  above  this  boundary  the  positive 
work  done  by  the  reaction  force  can  no  longer  be  expected,  and  further 
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Increase  in  the  limit-cycle  amplitude  Is  Impossible.  On  the  other  hand, 
when  the  positive  work  limit  appears  as  a  lower  boundary,  the  locus  can 
cross  the  limit,  and  once  the  crossing  occurs  the  amplitude  of  the  blade 
vibration  suddenly  drops  to  a  low  value  at  which  the  synchronized  shedding 
condition  Is  lost.  In  this  case,  even  if  the  positive  work  limit  exists 
above  the  stable-synchronization  limit,  the  Jump  from  the  strong  singing 
state  to  the  weak  one  should  appear  at  this  crossing. 

Thus  the  particular  feature  in  the  singing-frequency-versus-f low- 
velocity  relation  is  Interpreted  by  the  self-exciting  model.  Further, 
provided  that  the  mode  shape  of  blade  vibration  and  its  mechanical  damping 
are  given,  the  quantitative  derivation  of  the  f requency-versus-veloci ty 
diagram  Is  also  possible,  at  least  in  principle.  Before  entering  into 
these  quantitative  details,  however,  the  assumptions  and  simplifications 
adopted  in  the  mathematical  development  should  be  refined  and  reinforced. 
Investigation  of  most  of  them  Is  left  for  future  work,  but  some  discussion 
is  offered  below. 

The  approximations  given  for  the  mathematical  expressions  of  the 
Influence  functions  (x)  and  Sg(x)  do  not  seem  to  alter  the  nature 
of  the  problem,  at  least  qualitatively.  Their  essential  function  as  factors 
determining  the  vortex-shedding  frequency  comes  from  the  fact  that  (x) 
has  a  finite  value  starting  with  x  -  0  ,  while  (x)  Is  zero  at  x  >  0 
and  has  a  finite  incubation  interval  before  the  rapid  increase  with  x 
and  later  diminishing;  and  these  characteristics  of  (x)  and  S^(x) 
are  kept  in  the  simplified  expressions  of  Equation  (II). 

Another  speculative  assumption  was  made  concerning  the  concentrating 
process  of  the  vorticity.  To  treat  the  vortex-shedding  process  In  this 
way  instead  of  solving  the  Navier-Stokes  equation  may  be  allowable  when 
the  Reynolds  number  Is  not  so  low  and  the  transportation  of  the  vorticity 
by  flow  Is  dominant  compared  with  the  diffusion  of  the  vorticity  due  to 
viscosity.  For  such  a  case,  the  most  reliable  method  is  to  trace  the  path 
of  each  vortex  element  shed  from  the  separation  points,  in  the  manner 
utilized  by  Rosenhead  for  the  Helmholtz  instability  problem.  This  seems 
within  the  capacity  of  the  ordinary  computer.  In  the  present  work,  however, 
this  numerical  process  was  replaced  by  a  largely  simplified  mathematical 
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model,  In  order  to  Introduce  Into  the  hydrodynamic  problem  -  as  generally 
as  possible  —  the  concepts  developed  in  the  non-linear  oscillation  theory, 
and  to  clarify  the  present  version  of  the  singing  phenomenon  as  self-excited. 
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CONCLUSIONS 

A  mode)  for  the  propeller-singing  phenomenon  considered  as  a  self- 
excited  oscillation  was  presented  to  interpret  the  step  and  Jump  character¬ 
istics  In  the  s inging-f requency-versus-f low-veloci ty  diagram. 

The  singing  system  was  simulated  by  a  "closed  loop"  composed  of  a 
blade  as  a  mechanical  vibration  system  anJ  the  vort ices-shedding  mechanism 
which  Is  responsible  for  the  shedding  of  vortices  in  synchronization  with 
blade  vibration.  The  numerical  calculations  made  for  several  types  of 
blade-vibration  mode  shape  showed  that  the  criteria  for  the  stability  of 
the  synchronized  shedding  of  vort ices »  together  with  the  criterion  for 
the  phase  relation  favorable  to  positive  work  done  by  the  hydrodynamic 
reaction  force  of  the  shed  vortices  upon  the  blade  vibration,  can  interpret 
the  step  and  Jump  phenomenon. 

In  order  to  Introduce  the  concepts  of  the  non-linear  oscillation 
theory  into  the  hydrodynamic  problem  as  generally  as  possible,  large  sim¬ 
plifications  were  made  In  the  mathematical  expressions,  and  consequently 
the  quantitative  detailed  treatment  required  for  each  particular  case  Is 
left  for  future  work. 
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FIGURE  I. 


THE  RELATION  OF  SINGING  FREQUENCY  VERSUS  FLOW  VELOCITY 
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FIGURE  3b.  SCHEMATIC  FIGURE  FOR  THE  VORTEX-SHE DOING  PROCESS 
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FIGURE  4b.  INFLUENCE  FUNCTIONS:  Suu(x),  lx).  S«*  (x)  AND  S.,,  (x)  FOR  CIRCULAR 
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FIGURE  5b.  TRANSFORMED  PLANE 
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FIGURE  8.  COORDINATE  SYSTEM  FDR  FLAT-PLATE  BLADE 


FIGURE  9o.  ASSUMED  VIBRATION  MODE  OF  BLADE 


FIGURE  9b.  EXAMPLE  OF  AUXILIARY  DIAGRAM  FOR  A=0.6/27r,  b=2,  S' =  1.32  AND  0=1,2,  3  CORRESPONDING 
TO  S/S»l.20,  I.IO  AND  1.05  RESPECTIVELY 
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FIGURE  9c.  EXAMPLE  OF  AUXILIARY  DIAGRAM  FOR  A-O.S/Zw,  b  =  2.  S'*l.08  AND  0=1,2,  3  CORRESPONDING 
TO  S'/S  =0.98,  0.90  AND  0.86  RESPECTIVELY 
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FIGURE  lOc.  EXISTENCE  REGION  OF  SELF-EXCITED  SINGING 
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